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Abstract. In the context of Mass Varying Neutrinos(MaVaNs) we study a model in which a scalar 
field is coupled to more than one species of neutrinos with different masses. In general, adiabatic 
models of non-relativistic MaVaNs are heavily constrained by their stability towards the formation 
of neutrino nuggets. These constraints also apply to models with more than one neutrino species, 
and we find that using the lightest neutrino, which is still relativistic, as an explanation for dark 
energy does not work because of a feedback mechanism from the heavier neutrinos. 
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1. INTRODUCTION 



Precise observations of the cosmic microwave background [H], 0], the large scale struc- 
ture of galaxies [3], and distant type la supernovae [4] have led to a standard model 
of cosmology in which the energy density is dominated by dark energy with negative 
pressure, leading to an accelerated expansion of the universe. 

A proposal^ to explain dark energy is the so-called mass varying neutrino (MaVaN) 
model 



mm, 



14, 15 



in which a lig ht scalar field couples to neutrinos, see also js], %, [13, flU, 

M H H IM IM lai S IM IM IM l23, lai . 



In this paper we discuss the suggestion that the lightest neutrino which can be rela- 
tivistic today may be responsible for dark energy. We find that there is some evidence 
that the relativistic neutrino will feel an instability towards the formation of neutrino 
nuggets. 

In the next section we briefly review the formalism needed to study mass varying 
neutrinos, in section 3 we discuss MaVaNs with a relativistic neutrino, and in section 4 
we conclude. 



2. MASS VARYING NEUTRINOS 

In the MaVaN model [Isi |6i 0] we introduce a coupling between neutrinos and a light 
scalar field, and the coupled fluid then acts as dark energy. In this way, the neutrino 
mass niy is generated from the vacuum expectation value (VEV) of the scalar field. The 
effective potential is defined by 



(1) 



where V(j,{^)is the scalar field potential, a is the scale factor, pv (my ) , a) is the neutrino 
energy density, and Py{my{(^),a) is the neutrino pressure. 

The energy density and pressure of the scalar field are given by the usual expressions, 

p^[a) = ^f + V^{^) and P^[a) = -V^{^). (2) 

Defining w = /bE/PoE to be the equation of state of the coupled dark energy fluid, 
where /^E = + denotes its pressure and Pde = pv + P(p its energy density, the 
requirement of energy conservation gives, 

Pde + 3HPde{1+w)=0, (3) 

where H = ^ and dots to refer to the derivative with respect to conformal time. Com- 
bining with Eq. ([3]), one arrives at a modified Klein-Gordon equation describing the 
evolution of 

^+2H<i)+a^Vj, = -a^l5{py-3Py), (4) 

where primes denote derivatives with respect to (' = d/d^) and /3 = '^^°f^^ is the 
coupling between the scalar field and the neutrinos. 

It can be quite instructive to look at the behavior of MaVaN models in the case of 
non-relativistic neutrinos /V — 0, such that Eq. ([U) takes the form 

V=Py+V^ (5) 

Naturalness suggests we pick a scalar field mass(Curvature scale of the potential) to 
be much larger than the expansion rate of the Universe, 

V" = (/3' + /32) + y; ^ ml » H\ (6) 

In this case, the adiabatic solution to the Klein-Gordon Eq. ^ applies As a conse- 
quence, the scalar field will sit in the minimum of its effective potential V at all times 

^'=p;+v;=/3pv+v;=o (?) 

MaVaNs models can become unstable on sub-Hubble scales m^^ < a/k < H^^ in the 
non-relativistic regime of the neutrinos, where the perturbations 5pv evolve adiabati- 
cally. 



In Ref. [|18|l (see also [|3d. 13 IL I32L I33L \34\) it is shown that the equation of motion for 
the neutrino density contrast in the regime m^^ < a/k < H^^ can be written as 



^^CDM^CDM + ^b^b 



(8) 



where 

n — n \ \ JL- ^^^^ _ 



Geff = Gll + r—-^ and a = -^^. (9) 



Since neutrinos interact through gravity as well as through the force mediated by the 
scalar field, they feel an effective Newton's constant Geff as defined in Eq. The force 
depends upon the MaVaN model specific functions /3 and V(j, and takes values between 
GandG(l+2/32M2) on very large and small length scales, respectively. 

In certain cases of strong coupling neutrinos suffer an instability towards clumping 
in which case they stop behaving as dark energy In Ref. [18] a criterion for the 

stability was developed. (^'^ '\~~^^fv^^^~~p^ ^vSv < ^^cdm^cdm + ^b^b- This can be 

recast in a more convenient form , Cly < ^m, where we have neglected 

the effect of baryons compared to cold dark matter and we have assumed the density 
contrasts of roughly the same order. 

From the considerations above one can establish a list of criteria that MaVaN models 
need to fulfill. This was done in Ref. [29] where it was stressed that for single field 
MaVaN models that satisfy adiabaticity, the right amount dark energy today, correct 
neutrino mass as well as stability cannot be simultaneously fulfilled. This has previously 
been stated by Refs. [11] and [18]. 

Hence it has been suggested, in the context of multiple scalar field models, that 
neutrinos may be stable towards clustering if our effective potential has two minima: 
A false minimum in which our universe sits and a true minimum. The offset between the 
two minima is then interpreted as the dark energy density. The model is implemented in 
SUSY to avoid problems with small scalar field masses. Stability is ensured by letting the 
lightest relativistic neutrino be responsible for the dark energy [16t] . Below we analyze 
this suggestion. 



3. IMAVAN IVIODEL WITH A RELATIVISTIC NEUTRINO 

We assume the scalar field couples to all light neutrino species for naturalness reasons. 
In the case of three hierarchical neutrino masses, one would naively assume that, as a 
result of the coupling, the two heavier neutrinos would become unstable to clustering, 
whereas the lightest would remain stable. In the following we will argue that this is 
not possible, since a feedback from the growth of the heavy neutrino perturbations will 
cause the relativistic neutrino perturbation to grow as well. The following equation(see 
Ref. llJ) explains this 

5pv = ^1 q^dqda£fo{q)'i'+5(^P{py-3Py), (10) 

where = + niyO^ and /o is the unperturbed Fermi-Dirac distribution with a small 
perturbation ^. The equation explains the growth of the neutrino perturbation and 
applies to both the interaction between a relativistic neutrino and a scalar field as well as 
that of a non-relativistic neutrino and a scalar field. 

We consider a system consisting of two heavy neutrinos and one relativistic neutrino 
each interacting with the same scalar field. The following list of events will take place. 
In the beginning, the relativistic neutrino will not feel the presence of the heavier ones. 



However, the heavy neutrinos will feel the coupling which will drive their perturbations 
5pv (heavy) to large values. In this way, the heavy neutrinos start clumping. 

As was demonstrated in Ref. [18], the scalar field perturbation is effectively propor- 
tional to the neutrino perturbation for the interaction with heavy neutrinos. This means 
that since 5pv (heavy) is growing, 5^ will also grow accordingly. 

Regarding the relativistic neutrinos, the perturbations 5pv(rel) will grow as indicated 
in Eq. [TOl The second term in the equation consists of three important contributions, 
firstly 50 which is growing (this is the same 50 as listed above since we only have 
one scalar field). Secondly we have a coupling which we for simplicity assume to be 
constant (in reality this will be a growing quantity for most cases). And finally we 
have (pv — 3Pv) which is a suppression factor of the order m/E - this factor will act to 
delay the growth of 5pv(rel). However, since 50 will continue its growth, the inevitable 
conclusion is that 5pv(rel) will eventually start to grow. Hence, there exists a type of 
feedback mechanisms between the heavy and the relativistic neutrinos. One could of 
course argue that we are exactly living in a transition regime when 5pv(rel) has still not 
turned unstable. However, that would require serious fine-tuning. 

A graphical illustration of the example above is given in Fig. [IJ This is done in the 
framework of a model with a Coleman-Weinberg type scalar field potential similar to 
the one presented in Refs.LIJ and [[l8|l . 

V^=Vo\og{l+k^f) (11) 

In order to avoid possible pathological behaviour, we choose a mass term slightly 
different than the one in Ref. [18], namely one that does not become infinite when the 
VEV of the scalar field goes to zero. However, it still behaves as 1/0 for small 0. 

mv = -^A0 + y^A202 + m2, (12) 

which can be derived from solving the mass matrix (A0,m^;m^,O). 

What happens is shown in Fig. [H where we can see that for the higher redshifts the 
density contrast of heavier neutrinos behaves moderately as predicted by GR. The cdm 
term in eq. [8]sources the slow growth of these heavy neutrinos. As their masses increase 
more and more the coupling term slowly takes over and becomes the dominant term 
in eq. [81 Eventually this leads to the unstable growth of their density contrast. What 
happens next is that once the growth of the heavy neutrinos enter the quasi-linear regime, 
immediately the effect can be seen on the growth of the relativistic neutrino density 
contrast. This starts blowing up, and a short while later the system of equations we are 
solving effectively breaks down, which can be seen by the unnatural strong growth of 
the cdm density contrast. 

This gives us a hint that as a result of a feedback mechanism, the fast growing behavior 
of the heavy neutrino density contrast causes the relativistic neutrinos to start clumping 
as well. 

Hence the neutrino scalar- field fluid will start acting as a cold dark matter component 
(clustering neutrinos) and hence cannot be attributed to dark energy. 

Now, one can of course argue that the rise in the density contrast of the relativistic 
neutrino species only happens once the heavier ones have turned non-linear. In this 
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FIGURE 1. Density contrasts plotted as a function of redshift for a system consisting of one light and 
two heavy MaVaN neutrinos each interacting with the same scalar field. The scale is k = 0.1 Mpc^' and 
we choose the current neutrino masses my (rel) = 10^^ eV and my (heavy) = 0.3 eV (Note that the choice 
of current neutrino masses does not affect the result qualitatively). The solid line is cdm-density contrast, 
the dotted line is the light neutrino density contrast, and the dashed line is the heavy neutrino density 
contrast. The heavy neutrinos grow essentially as cdm until the coupling becomes large enough for the 
instabilities to set in. The light neutrino is still relativistic, and its density contrast oscillates as acoustic 
waves. However, due to a feedback mechanism, the relativistic neutrino density contrast tracks that of the 
non-relativistic neutrino around the time the growth of the heavy neutrino perturbations become quasi- 
linear. I.e. both neutrino species will clump. Note that the cdm perturbations also blow up at late times. 
This is an effect of the system of differential equations breaking down as all parameters go to infinity. 

regime the linear code does not apply. However, from a close-up look at the data, we 
emphasize that the rise of the density contrast happens in the quasi-linear regime of the 
heavy neutrino perturbations, where the code does still apply. 

The reason that the relativistic neutrino is able to clump is that it will acquire an 
effective mass, thus it cannot be regarded as a relativistic particle. Unfortunately, we 
cannot use conventional bounds to constrain this effect for the following reason: Once 
the evolution of the non-relativistic neutrinos becomes non-linear, the whole system of 
equations we are solving, starting with the modified KG equation breaks down. This has 
the effect that all current bounds are no longer valid, as these are established in the linear 
regime. 

4. CONCLUSION 

Single scalar field models can be used to explain late-time acceleration in the MaVaN 
scenario. However, in general using these potentials leads to instabilities towards neu- 
trino bound states unless certain criteria are relaxed. 

Accordingly it has been suggested to include an extra scalar field in the treatment. This 
has some very nice features and is easily capable of obtaining late-time acceleration as 
well as Q-oE = OJ today. However, one drawback is the need for the lightest neutrino to 




be relativistic today. As was explained above, the feedback mechanism will eventually 
cause the relativistic neutrino to start clustering and hence the coupled fluid will cease 
to act as dark energy. 
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